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ORIGINAL ARTICLE 

Assessment of Maraviroc Exposure-Response 
Relationship at 48 Weeks in Treatment-Experienced 
HIV-1 -Infected Patients in the MOTIVATE Studies 

P Jacqmin 1 , JR Wade 1 , B Weatherley 2 , E Snoeck 1 , S Marshall 2 and L McFadyen 2 

Efficacy exposure-response relationships of the CCR5 antagonist maraviroc were evaluated across two phase III clinical trials. 
This post-hoc analysis used 48-week efficacy data from 841 treatment-experienced patients infected with CCR5-tropic human 
immunodeficiency virus type 1 (HIV-1), identified by the enhanced sensitivity Trofile assay. Probability of treatment success 
(viral RNA <50 copies/ml) was modeled using generalized additive logistic regression, testing exposure, clinical, and virologic 
variables. Prognostic factors for treatment success (in decreasing order of Akaike information criterion (AIC) change) were: 
maraviroc treatment, high-weighted overall susceptibility to background treatment, absence of an undetectable maraviroc 
concentration, high baseline CD4 count (BCD4), low viral load (VL), race (other than black), absence of non-R5 baseline tropism 
(BTRP), and absence of fosamprenavir (FPV). No concentration-response relationship was found with treatment (maraviroc 
vs. placebo) and presence/absence of undetectable maraviroc concentration (adherence marker) in the model. The maraviroc 
doses studied (300 or 150mg with potent CYP3A4 inhibitors once (q.d.)Awice daily (b.i.d.)) deliver concentrations near the top 
of the concentration-response curve. 
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Maraviroc is a reversible and selective antagonist of the 
human chemokine receptor CCR5. 1 It is the only entry inhibi- 
tor that has been approved in combination with other anti- 
retroviral agents to treat patients infected with CCR5-tropic 
human immunodeficiency virus type 1 (HIV-1). This approval 
was based on two studies, MOTIVATE 1 (NCT00098306) and 
MOTIVATE 2 (NCT00098722). 2 The studies were conducted 
in heavily treatment-experienced patients who not only had 
complex medical and treatment histories but also were tak- 
ing many concomitant medications (including optimized 
background therapy) with potential pharmacokinetic (PK) 
and pharmacodynamic (PD) interactions, which could affect 
exposure-response relationship and clinical outcome. 

The term "curse of dimensionality" was coined by Richard 
Bellman to describe the problem caused by the exponential 
increase in volume associated with adding dimensions to a 
mathematical space. 3 In the context of building a model to 
predict clinical response to antiretrovirals and other anti-infec- 
tives, including maraviroc, the curse of dimensionality means 
that as the number of potential predictors increases it becomes 
harder to find the best model. By convention, established PK/ 
PD analysis methods/models assume that concentration is 
the main driver of response. This might not necessarily be the 
case, especially if the dose(s) deliver concentrations toward 
the top (or bottom) of the concentration-response curve. It 
has long been recognized that clinical response is not solely 
dependent on PK, but there are many other factors relating 
to the status of the HIV-positive patient that also play a role. 45 
This has clearly been shown to be the case for maraviroc par- 
ticularly when focusing on overall virologic success. 2 For this 
reason, generalized additive models (GAMs) were employed 



to characterize the influence of prognostic factors (including 
exposure parameters) on sustained virologic response and 
incidence of anemia in patients with chronic hepatitis C 6 and 
in HIV-1 -infected patients in the etravirine phase III clinical 
studies. 7 GAMs allow potential predictors to enter linearly or 
nonlinearly, as appropriate for each case. 

In the MOTIVATE 1 and MOTIVATE 2 studies, after 48 
weeks' treatment, more patients receiving maraviroc once 
(q.d.) or twice daily (b.i.d.) with optimized background treat- 
ment (OBT) had HIV-1 RNA levels <50 copies/ml than those 
treated with OBT alone (MOTIVATE 1 : 42 and 47% vs. 16% 
MOTIVATE 2: 45 and 45% vs. 18%, respectively). 2 Sub- 
group analyses of the data pooled from both the MOTIVATE 
studies were performed and the results of the multivariate 
logistic regression modeling showed that virologic response 
at week 48 of maraviroc administration was significantly 
related to race, viral load (VL) at screening, CD4 cell count, 
and OBT containing enfuvirtide (first use). 8 Subgroup anal- 
ysis of the MOTIVATE studies was extended to examine 
weighted OBT susceptibility score category (WOBTSSC), 
derived by combining genotypic or phenotypic resistance 
data with prior drug use, to assess OBT activity. 9 This logis- 
tic regression analysis showed that genotypic and pheno- 
typic weighted scores were better at predicting response 
than counting active drugs. However, no measures of mara- 
viroc PK/exposure were included in these published sub- 
group analyses. 

A previously reported prespecified GAM analysis (per- 
formed at 24 weeks) for the endpoint of probability of fail- 
ure (HIV RNA >50 copies/ml; missing = failure) found a 
sigmoid-type exposure-response relationship in addition to 
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other disease and virologic factors. 10 This analysis did not, 
as planned, use pill counts as an adherence measure in the 
models because available pill-count data were thought to 
be unreliable. The objectives of this post-hoc 48-week GAM 
analysis were to assess exposure-response relationship and 
other possible predictive factors, as before. However, the 
concentration-based variables were derived without PK sam- 
ples below the limit of quantification (BLQ) and an additional 
covariate, presence/absence of one or more such samples 
(assumed to result from >3 consecutive/clustered missed 
doses) was included to account for poor adherence. 

RESULTS 

The analysis was performed using an endpoint where all 
patients with missing 48-week efficacy data (i.e., missing or 
discontinued to week 48) were treated as failures (MD = F) 
(total n = 841 including 178 placebo-treated patients). Informa- 
tion on prognostic factors by treatment group for this dataset 
is shown in Table 1 . Patients with one or more missing cat- 
egorical or continuous prognostic factors listed in Table 1 were 
removed before running the step-wise GAM process, which 
left 782 patients. Twenty-five maravi roc-treated patients had no 
usable maraviroc concentration data and therefore were also 
excluded. After selection of the final model, those patients with 
missing data for covariates that had not been selected for inclu- 
sion in the model were returned to the final dataset (n = 801) 
before re-running the final model. 

Exploratory exposure-response analysis 

A graphical analysis of the relationship between average 
plasma concentration (C ), minimum plasma concentration 
(C min ), and effective constant concentration (ECC = the con- 
stant concentration giving the same average viral inhibition 
as the time varying concentration) demonstrated, not unex- 
pectedly, that they were highly correlated (data not shown). 

A graphical analysis of the virologic success (VL <50 cop- 
ies/ml) rate by ECC exposure quartile indicated a marked 
increase between the placebo rate (1 5%) and the first quartile 
(39%), a less pronounced increase between the first and the 
second quartiles (50%), and a slight further increase from the 
second to the fourth quartile (57%) (Figure 1). This appar- 
ent relationship between the success rate and exposure was 
further investigated by logistic regression with GAM to allow 
adjustment for other factors in making this assessment. 

GAM analysis of virologic success 
Supplementary Table S1 online shows the prognostic fac- 
tors (and relevant coding) for testing in the GAM analysis. 

The results of the automated step-wise searches for the prog- 
nostic factors of virologic success (VL <50 copies/ml; MD = F) 
at week 48 are given in Table 2. The initial search (Model 1) 
tested the following factors: baseline CD4 count (BCD4); base- 
line VL (BVL); treatment group (TRT); baseline tropism (BTRP); 
WOBTSSC; protease inhibitor group (PIGRP); presence of rito- 
navir, saquinavir, tipranavir, atazanavir, indinavir, lopinavir, and 
fosamprenavir (FPV) in the OBT; one or more concentrations 
BLQ; age; sex; and race. The retained prognostic factors were: 
BCD4, BVL, TRT, WOBTSSC, RACE, and BLQ. 



Table 1 Description of the prognostic factors used in the analysis by 



treatment group 



Parameters 


Placebo 


q.d. 


b.i.d. 


Number of patients 


178 


335 


328 


Mean age (years) 


45.7 


45.8 


46.6 


ft /I 1 _ /o/ \ 

Male (%) 


87.6 


86.9 


88.7 


vvnite \/o) 


84.8 


80.3 


85.1 


Black (%) 


12.4 


17.9 


12.2 


Asian (%) 


0.6 


0.6 


0.9 


Other (%) 


2.2 


0.9 


1 .8 


Success (VL <50 copies/ml) (%) 


15.2 


44.5 


47.9 


Failure (VL >50 copies/ml) (%) 


84.8 


55.5 


52.1 


Virologic failure (%) 


74.7 


46.6 


45.1 


Nonvirologic failure (%) 


1 0.1 


9.0 


7.0 


Median average concentration (ng/ml) 


0.0 


108.0 


205.7 


Median minimum concentration (ng/ml) 


0.0 


32.7 


96.0 


Median effective constant concentration (ng/ml) 


0.0 


65.1 


171.2 


Patients with BLQ (%) 


0.0 


7.2 


3.7 


Patients with no exposure information (%) 


0.0 


3.9 


3.7 


Median baseline CD4 count (cells/mm 3 ) 


178.0 


173.5 


180.3 


Mean log 10 (baseline VL) (copies/ml) 


4.9 


4.9 


4.9 


Patients with R5 tropism at baseline (%) 


92.1 


93.7 


95.1 


Patients with X4 or dual-mixed tropism at 


6.2 


3.6 


3.0 


baseline (%) 








Patients with non-readable tropism at baseline (%) 


1 .7 


2.7 


1 .8 


Patients with baseline WOBTSS, 0-0.5 (%) 


36.0 


35.2 


32.0 


Patients with baseline WOBTSS, 1-1 .5 ( %) 


33.7 


34.6 


35.1 


Patients with baseline WOBTSS, 2-2.5 (%) 


25.3 


24.8 


26.2 


Patients with baseline WOBTSS, >3 (%) 


5.1 


5.4 


6.7 


Patients treated with ritonavir (%) 


86.0 


85.7 


85.7 


Patients treated with saquinavir (%) 


1 1 .2 


15.2 


10.7 


Patients treated with tipranavir (%) 


12.9 


14.0 


14.9 


Patients treated with atazanavir (%) 


15.2 


16.7 


14.6 


Patients treated with indinavir (%) 


1.7 


2.7 


2.4 


Patients treated with lopinavir (%) 


25.3 


30.4 


35.7 


Patients treated with fosamprenavir (%) 


30.9 


20.0 


23.2 


Patients treated with bPI (%) 


47.8 


55.8 


53.7 


Patients treated with bFPV (%) 


24.2 


14.9 


17.1 


Patients with treatments other than bPI and 


28.1 


29.3 


29.3 



bFPV (%) 

bFPV, boosted fosamprenavir; b.i.d., twice daily; BLQ, below the limit of 
quantification; bPI, boosted protease inhibitor; q.d., once daily; VL, viral 
load; WOBTSS, weighted optimized background therapy susceptibility score 
category. 

Model 2 included all the previous prognostic factors except 
for TRT, but included ECC instead. In this search, the same 
prognostic factors as in Model 1 were retained with the addi- 
tion of ECC (probably in place of TRT), BTRP, and PIGRR 
However, the Akaike information criterion (AIC) was higher 
than in Search 1 (809.015 vs. 802.506). 

Model 3 tested all the previous prognostic factors including 
TRT and ECC with identical results to the first search indicat- 
ing that TRT was more predictive than ECC. 
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35.4 


50 
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47 
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47 
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332.0 
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56 


294.5 


57 
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Figure 1 (a,b,c) 
concentration; C 



Graphical and (d) tabular display of quartile analysis of response rates by C avgJ C min , and ECC. 
, minimum plasma concentration; ECC, effective constant concentration; VL, viral load. 



C avg , average plasma 



In Model 4, WOBTSSC was removed from the list of 
prognostic factors tested. The resulting model indicated that 
BTRP and FPV only partially accounted for response differ- 
ences explained by WOBTSSC, as AIC was markedly higher 
(860.974). 

Model 5 investigated the use of OBT overall susceptibility 
score category (OSS) instead of WOBTSSC. OSS and BTRP 
were selected but AIC (824.331) was higher than when 
WOBTSSC was included. 

Model 6 included all previous prognostic factors in the 
search with ECC as the concentration variable. The result 
was identical to Models 1 and 3, confirming that TRT and 
WOBTSSC are better prognostic factors than ECC and OSS, 
respectively. 

TRT may be a better prognostic factor than ECC because 
the information about placebo and/or the maraviroc regimen 
(q.d. and b.i.d.) contained in TRT was similar or superior 
to ECC in explaining study outcome. This raised the ques- 
tion of whether the distinction between maraviroc q.d. and 
b.i.d. is informative or necessary. Therefore prognostic factor 
TRT2, which contains only placebo and maraviroc (q.d. and 
b.i.d. combined) categories, was added to the search. In this 
search (Model 7), TRT2 was selected instead of TRT and AIC 
was slightly lower (800.540 vs. 802.506) than when TRT was 
selected, suggesting no difference between maraviroc q.d. 
and b.i.d. when the other prognostic factors were taken into 
account. 

In the current analysis, BTRP was not consistently selected, 
whereas it was a significant prognostic factor in all models 
tested in the previous MOTIVATE 24-week analyses 10 where 
subjects were screened with the original Trofile assay. This is 



probably because the signal is weaker in this analysis with 
fewer patients with non-R5 virus at baseline as some were 
screened out by the more discriminating enhanced sensitiv- 
ity Trofile assay. Therefore, prognostic factor BTRP2, catego- 
rizing only maravi roc-treated patients for BTRP was derived 
and added in Search 8. BTRP2 and FPV were selected in 
addition to the factors already selected in Search 7, but with 
a slightly lower AIC (798.91 7 vs. 800.540) and with very small 
contributions to the AIC changes (Table 2; Figure 2). This is 
consistent with a small differential effect of BTRP in predict- 
ing outcome in subjects in maraviroc groups as compared 
with placebo. 

As ECC was not selected in the previous search, Searches 

9 and 10 were performed to check whether C min or C avg may 
be better predictors of response than ECC. The results again 
showed that neither C . nor C was more informative than 

min avg 

TRT2 and BLQ. 
In the best model obtained at this stage (Searches 8, 9, and 

10 yielding identical models), BLQ was the only concentration- 
related parameter retained as a prognostic factor. In order to 
test whether BLQ was preventing other exposure variables 
from entering the model, BLQ was removed and two new 
searches were conducted: one with all prognostic factors and 
TRT (Search 11) and another with all prognostic factors and 
TRT2 (Search 12). ECC was selected in the search with TRT, 
but not in the search with TRT2. The AlCs were higher than 
that obtained in Search 8 (835.729 and 835.310 vs. 798.917). 
This indicated that BLQ was the most informative concentra- 
tion-related parameter. It is important to note that the concen- 
tration variables were derived in the PK analysis with BLQ 
excluded from the PK dataset, thus exposure variables were 
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Table 2 Generalized additive models and Akaike criteria obtained at the end of automated step-wise searches for prognostic factors of virologic success 
(VL <50 copies/ml) at week 48 in HIV patients from MOTIVATE 1 and MOTIVATE 2 studies (final Model 8 in bold) 



Search 


Test 


Model 




AIC 


1 


TRT in absence of exposure 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + TRT + WOBTSCC + RACE + BLQ 


802.506 


2 


ECC in absence of TRT 


I ill //"*i 1 1~ r\ 

HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + lo(ECC) + BTRP + WOBTSCC + PIGRP + RACE + BLQ 


809.015 


3 


ECC in presence of TRT 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + TRT + WOBTSCC + RACE + BLQ 


802.506 


4 


Absence of WOBTSSC 


i ill 1 1~ r\ 

HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + TRT + BTRP + FPV + RACE + BLQ 


860.974 


5 


OSS instead of WOBTSSC 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + TRT + BTRP + OSS + RACE + BLQ 


824.331 


6 


ALL(ECC) 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + TRT + WOBTSCC + RACE + BLQ 


802.506 


7 


ALL(ECC) + TRT2 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + WOBTSCC + RACE + BLQ + TRT2 


800.540 


8 


ALL(ECC) +TRT2 + BTRP2 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + WOBTSCC + FPV + RACE + BLQ + BTRP2 +TRT2 


798.917 


9 


ALL(CAVE) + TRT2 + BTRP2 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + WOBTSCC + FPV + RACE + BLQ + BTRP2 + TRT2 


798.917 


10 


ALL(CMIN) +TRT2 + BTRP2 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + WOBTSCC + FPV + RACE + BLQ + BTRP2 + TRT2 


798.917 


11 


ALL(ECC) + BTRPT2 - BLQ 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + TRT + s(ECC, 3) + WOBTSCC + RIT + TIP + RACE + BTRP2 


835.729 


12 


ALL(ECC) + TRT2 + BTRP2 - BLQ 


HIVCV50 ~ 


log(BCD4) + lo(log(BVL)) + WOBTSCC + FPV + RACE + BTRP2 + TRT2 


835.310 



ALL = all factors listed as additional prognostic factors except OSS. + (plus): Adds an additional named factor to those provided for the search. - (minus): 
Removes named factor from those provided for the search. 

AIC, Akaike information criterion; BCD4, baseline CD4 count; BLQ, below the limit of quantification; BTRP, tropism at baseline; BVL, baseline viral load; ECC, 
effective constant concentration; FPV, fosamprenavir; HIV, human immunodeficiency virus; OSS, overall susceptibility score; PIGRP, protease inhibitor group; 
RIT, ritonavir; TIP, tipranavir; TRT, treatment group; VL, viral load; WOBTSSC, weighted optimized background treatment susceptibility score category. 



complementary to BLQ. Search 8 was therefore chosen as 
the final model. The sequence of prognostic factors entering 
this model during the search is given in Figure 2. In Table 3, 
the contribution of the individual prognostic factors to the AIC 
change derived by removing the parameters one at a time is 
shown in descending order of AIC change together with the 
direction of change. In order of importance, the prognostic fac- 
tors associated with success were: treatment with maraviroc 
(vs. placebo), higher WOBTSSC, absence of a maraviroc BLQ 
(good adherence), higher CD4 count, lower BVL, race (other 
than black) with only very minor changes in AIC associated 
with tropism, and absence of FPV. Odds ratios for categorical 
and linear terms are included to indicate the magnitude of their 
effects. 

Figure 3 depicts the uncertainty of the predicted percentage 
of success in the investigated patient population. The narrow 
distributions indicate that the model parameters were estimated 
with acceptable uncertainty. The predictive performance of the 
GAM to simulate trial outcomes was also successful as evalu- 
ated by a visual predictive check and where the trial simula- 
tions accounted for both the parameter uncertainty (from 1 ,000 
bootstraps) as well as the residual error (data not shown). 

Figure 4 represents the predicted probability of success 
(VL <50 copies/ml) at week 48 using the final GAM (Model 
8 in Table 3) showing the effects of the exposure covariates. 
Figure 4 (top row, left) shows mean predicted (and observed) 
probability of success at week 48, with 1 5% probability of suc- 
cess for placebo vs. 46% for maraviroc q.d. and 49% for mara- 
viroc b.i.d. In the two-treatment category model (top, right), the 
predictions (and observed values) are 47.4% for maraviroc 
(q.d. + b.i.d.) vs. 15% for placebo. In the bottom row (right), tak- 
ing the BLQ variable into account for the maraviroc group (q.d. 
+ b.i.d.), the predicted probability of success increased to 50% 
for those without an observed BLQ whereas the probability of 
success for those with one or more BLQ values on maraviroc 
was only 3%, lower than that for placebo at 1 5% (bottom, left). 



BCD4 
TRT2 

-} WOBTSSC 
-> BLQ 
-> log(BVL) 
-> RACE 

BCD4 log(BCD4) 
log(BVL) s(log(BVL)), 2 

s(log(BVL, 2)) s(log(BVL)), 3 

s(log(BVL, 3)) -fr lo(log(BVL)) 
-» BTRP2 
FPV 



800 850 900 950 1,000 1,050 
Akaike information criterion 

Figure 2 Results of the automated step-wise search for prognostic 
factors of virologic success (viral load (VL) <50 copies/ml) at week 
48 in human immunodeficiency virus patients from the MOTIVATE 
studies (Search 8 in Table 2). BCD4, baseline CD4 count; BLQ, 
below the limit of quantification; BVL, baseline VL; TRT, treatment 
group; WOBTSS, weighted optimized background therapy 
susceptibility score category. a, indicates addition of variable a; 
a^b, indicates the replacement of variable a with b. BTRP, baseline 
tropism; FPV, fosamprenavir. 

DISCUSSION 

The many potential predictors and the resulting curse of 
dimensionality present in the treatment of HIV infection are 
challenging in terms of both the range of potential predictors 
and the choice of methodology for finding the best model that 
includes predictors of clinical success. Different analysis meth- 
ods are available when analyzing data with multiple potential 
predictors, these include tree analysis, 11 neural networks, 12 
and the more traditional approach of regression modeling. 13 
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Table 3 Contribution of each prognostic factor to the AlC in the final model (from Search 8 in Table 2) obtained by removing factors from Model one at a time 



Test (direction of relationship to success) 


Odds ratio 


AlC 


AAIC 


Reference 




811.383 


NA 


TRT2 (t maraviroc vs. placebo) 


8.27; maraviroc vs. placebo 


893.855 


82.472 


WOBTSSC (T with higher WOBTSSC) 


2.35, 5.84, 9.62; WOBTSSC 1 , 2, 3 vs. 0 


872.886 


61.503 


BLQ (iwith BLQ) 


0.019; BLQ vs. non-BLQ 


851 .089 


39.706 


BCD4 (T with higher log CD4) 


3.75; for tenfold increase 


841.691 


30.308 


BVL (i with higher log viral load, nonlinear) 




834.135 


22.752 


RACE (i in black patients) 


0.35; black vs. white 


821.661 


10.278 


BTRP2 (i with non-R5 in maraviroc group) 


0.35; non-R5 vs. R5 


813.053 


1.670 


FPV (i in patients on fosamprenavir) 


0.74; FPV present vs. FPV absent 


811.435 


0.052 



Higher AAIC indicates higher contribution. The direction of the relationship is also shown together with odds ratios. 

AlC, Akaike information criterion; BCD4, baseline CD4 count; BLQ, below the limit of quantification; BTRP, tropism at baseline; BVL, baseline viral load; FPV, 
fosamprenavir; NA, not appropriate; TRT, treatment group; WOBTSSC, weighted optimized background treatment susceptibility score category. 



In the present analysis, logistic regression analysis with GAMs 
was the method chosen to analyze the VL data as it has been 
successfully used previously in the analysis of treatment suc- 
cess in chronic hepatitis C-infected patients and in HIV. 7 ' 810 

The current analysis included PK measures as poten- 
tial predictors but cannot be considered as a classical PK/ 
PD analysis. PK/PD analysis of VL endpoints in treatment- 
experienced HIV-1 -infected patients with OBT is complex, 
mainly because (i) the investigated drug is given in combina- 
tion with various OBT combinations for which the PK and PD 
interactions are not fully understood; (ii) left censoring of the 
HIV-1 RNA levels due to the lower limit of quantification; (iii) 
numerous strong nonlinear covariate effects related to the 
disease state; and (iv) dropout. In the present analysis, VL 
was transformed into a binary endpoint (success or failure) 
and logistic GAMs were used to evaluate the effect of maravi- 
roc exposure and possible interacting prognostic factors. The 
binary endpoint used for the GAM approach partly resolves 
the left censoring of the VL measurements and allows the 
combination of prognostic factors as categorical, linear, or 
nonlinear terms without making too many assumptions about 
the shape of the factor-effect relationship. With respect to the 
binary endpoint, all patients who discontinued for any reason 
before 48 weeks were regarded as treatment failures thereby 
dealing with dropout to some extent. A sensitivity analysis 
without nonvirologic failure patients led to the selection of 
identical prognostic factors (data not shown). 

The results of this analysis of maraviroc at the 48-week 
milestone indicate, in addition to the well-recognized prognos- 
tic factors such as BCD4-positive T-cell count, BVL, baseline 
WOBTSS, other factors including treatment vs. placebo group, 
presence of a maraviroc concentration BLQ in maraviroc- 
treated patients, and race (Caucasians and others vs. black 
race) were prognostic factors of virologic success. In addition, 
virus BTRP in maravi roc-treated patients and absence of FPV 
in the OBT were also identified as prognostic factors, but with 
a very small AlC change. These results are generally in agree- 
ment with other analyses of the MOTIVATE data. 28-10 

The finding that a binary treatment variable (placebo vs. 
maraviroc and no difference between maraviroc q.d. and 
b.i.d.) was preferred to concentration variables, suggests that 
success is similar across the concentration range covered 
by both the q.d. and b.i.d. treatments and dose levels in the 
present data set after accounting for other prognostic factors. 
It should be noted that -70% of the patients were comedi- 
cated with boosted Pis, which increase maraviroc area under 
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Figure 3 Visual predictive check of the generalized additive model 
for virologic success (viral load <50 copies/ml at week 48) in the 
MOTIVATE studies by dose group. The central red line represents the 
observed success rate in the studies. The side orange lines represent 
the 95% confidence intervals of the predictions, b.i.d., twice daily; FAS, 
full analysis set; OBT, optimized background therapy; q.d., once daily. 
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Figure 4 Predicted probability of success (VL <50 copies/ml) at week 48 using the final generalized additive model (Model 8 in Table 2) in 
patients from the MOTIVATE studies. Top row is by TRT group and TRT2 group, bottom row is by BLQ category, b.i.d., twice daily; BLQ, below 
the limit of quantification; MVC, maraviroc; q.d., once daily; TRT, treatment group: placebo, MVC q.d., MVC b.i.d.; TRT2, treatment group: 
placebo, MVC q.d. and b.i.d. combined; VL, viral load. 



the concentration vs. time curve (AUC) by two- to eightfold 
(dependent on the PI) because of CYP3A4/P-glycoprotein 
inhibition. Maraviroc doses used in the MOTIVATE studies 
appear to have delivered concentrations high on the E max 
(inhibitory maximum effect) exposure-response relation- 
ship consistent with maraviroc monotherapy dose-ranging 
studies 14 and the modeling of these data. 15 

In contrast to the previous analysis of the MOTIVATE 
24-week efficacy data, 10 none of the concentration-based 
exposure variables (C avg , C min , and ECC) were significant fac- 
tors in this analysis, in the presence of the treatment group 
and maraviroc BLQ. The earlier identification of a sigmoid rela- 
tionship for concentration-based exposure variables 10 is highly 
dependent on the presence of the placebo group (zero mara- 
viroc concentration). This can also be seen from the quartile 
analysis (Figure 1) in the present study, with fairly flat relation- 
ships across the wide maraviroc exposure ranges but a sharp 
increase in success from the placebo group to the first quartile. 

Another important learning is that inclusion of BLQ values 
(set to 0.5 x lower limit of quantification) derived from sparse 
PK sampling may introduce a bias in the population PK mod- 
eling if adherence is not independently identified and utilized 
in the analysis. This may result in the derivation of incorrect 



and/or misleading exposure-response relationships. Having 
a maraviroc BLQ associated with the dosing regimens used 
in the MOTIVATE studies is deemed a good surrogate for a 
treatment interruption of at least 2-3 days. This is based on 
maraviroc PK from phase l/l la studies where concentrations 
were measurable 48-72 h after oral doses of 300 mg q.d. 
or 100mg b.i.d. and the slow terminal elimination phase of 
maraviroc seen after intravenous dosing with 30 mg (terminal 
half-life of 13h). 16 In this analysis, we have found using BLQ 
as a separate prognostic surrogate factor for adherence and 
not including BLQ in the PK analysis provides a clear under- 
standing of the most likely prognostic factors. Our analysis 
indicates that cases of maraviroc BLQ identify poor adher- 
ence (probably to all antiretrovirals) as the response rate is 
only 3 vs. 1 5% for the placebo treatment group. 

Consecutive/clustered missed doses are the most likely pat- 
tern of poor adherence to cause treatment failure/viremia in 
HIV-1 -infected patients. Only 2 of the 36 patients with >1 BLQ 
values in the data set had a VL <50 copies/ml at week 48 with 
only 4 having discontinued for reasons unrelated to virologic 
failure. The observed imbalance in the number of patients with 
BLQ observations in the maraviroc q.d. (7.2%) vs. b.i.d. 
(3.7%) arms may be attributed to one or more of the following: 
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(i) patients in the q.d. arm received half the daily maraviroc 
dose of the b.i.d. arm (maraviroc 1 50/300 mg depending on 
OBT vs. 300/600 mg); (ii) the consequences of one missed q.d. 
dose is equivalent to two to three missed b.i.d. doses; 17 (iii) the 
active dose in the double-dummy dosing scheme for q.d. was 
given in the evening and there is some evidence to suggest 
that patients prefer to take medication in the morning. 18 

PK measurements were only performed up to and including 
the 24-week visit and this may be regarded as a shortcoming 
of the present analysis. However, the first 12-24 weeks of 
treatment are the most important in reducing viral RNA and 
when adequate treatment and very high adherence to treat- 
ment are critical. This is supported by the similar response 
rates seen in the maraviroc arms at 24- and 48-week time- 
points for <50 copies/ml viral RNA. 2 Therefore, the informa- 
tion and adherence behavior collected early in the study 
can probably be safely extrapolated to later timepoints for 
patients remaining in the study from 24 to 48 weeks. The big- 
gest drawback of the PK/PD analysis in these studies is that 
the accuracy of the PK/exposure covariates relies entirely on 
the accuracy of the patient-reported dosing information gath- 
ered at each visit. 19 Longitudinal single or even duplicate PK 
measurements at multiple clinic visits dependent on patient- 
reported dosing history are not an efficient way of determin- 
ing exposure-response relationships nor do they indicate 
the need for concentration monitoring. Electronic monitoring 
is a much more informative way of capturing adherence on 
an ongoing basis including between clinic visits, especially 
in HIV where sustained high adherence is required. 19 Pill- 
count data are probably the least informative way to measure 
adherence, especially where clustered missed dosing may 
be critical, as in HIV therapy. Attempts were made to use pill- 
count data in the 24-week analysis of the MOTIVATE studies. 
Pill counts were, however, not pursued in the step-wise GAM 
because a univariate analysis showed a bell-shaped relation- 
ship with the highest efficacy centered on 100%. Pill counts 
>100% arise when patients do not bring back any pills. In this 
scenario, it is often assumed that patients have good adher- 
ence while indeed this could be indicative of poor adherence. 

FPV in the OBT is somewhat over-represented in the 
placebo group (Table 1) and this may be responsible for its 
appearance as a minor covariate. It has been shown across 
a number of studies of anti retroviral drugs in treatment-naive 
patients that black patients generally show a worse response 
than other races. A recent meta-analysis across six studies 
for four antiretrovirals showed a 48% lower chance of achiev- 
ing VL <50 copies/ml at 96 weeks in black patients. 20 The 
researchers questioned whether this is because of adher- 
ence issues or other factors, such as differences in PK or 
psychosocial barriers. In the MOTIVATE PK covariate analy- 
sis, black patients appeared to have slightly higher PK con- 
centrations than other patients. 21 

In conclusion, provided datasets are large enough, the 
use of GAM-type modeling is an efficient method of assess- 
ing multiple covariates with possible linear and nonlinear 
relationships (including exposure variables) in landmark 
analyses. However, it is also important to carefully consider 
how adherence is factored into population PK with sparse 
sampling and how this is subsequently used in exposure- 
response assessments. 



It appears that the maraviroc doses used in the MOTIVATE 
studies delivered concentrations high on the E max exposure- 
response relationship. Therefore, no concentration-based 
exposure-response relationships were identified once other 
prognostic factors, including treatment group (maraviroc vs. 
placebo) and adherence (presence/absence of measurable 
maraviroc concentrations), were taken into account. 

METHODS 

Studies analyzed. This analysis used 48-week efficacy data 
from two previously reported 2 placebo-controlled phase lib/ 
III studies (MOTIVATE 1 and MOTIVATE 2) of maraviroc and 
OBT in treatment-experienced patients infected with CCR5- 
tropic HIV-1. The study protocols were approved by institu- 
tional review boards or independent ethics committees at 
each study center. Patients reported as being infected with 
CCR5-tropic HIV-1 in the MOTIVATE 1 and 2 studies were 
rescreened with the enhanced sensitivity Trofile assay 22 (total 
n = 841 with 208 censored from the original 1 ,049 for being 
non-CCR5). Maraviroc was given at a dose of 150mg q.d. 
or b.i.d. when administered with Pis (except tipranavi ^rito- 
navir) or delavirdine, otherwise 300 mg (q.d. or b.i.d.) was 
administered. 

VL was measured using the Amplicor HIV-1 monitor v1.5 
assay (Roche Diagnostics, Basel, Switzerland). OSS to the 
selected OBT (except for enfuvirtide) was assessed by geno- 
typic and phenotypic sensitivity (Phenosense GT assay; 
Monogram Biosciences, San Francisco, CA) with a simple 
binary score of 1 for full sensitivity and 0 for any reduced sus- 
ceptibility. For enfuvirtide, gp41 sequencing was performed 
(British Columbia Centre for Excellence in HIV/Aids, Vancou- 
ver, British Columbia, Canada). The WOBTSS (genotypic) 
excluded agents recycled from the pre-study regimen. All 
agents to which the screening virus was ranked as "sensi- 
tive" scored 1 .0 except for nucleoside reverse transcriptase 
inhibitors, which scored 0.5. Pis with an "intermediate" rank- 
ing (whether boosted or not) scored 0.5, whereas all other 
drugs ranked as "intermediate" scored 0. 

Exposure variables. Maraviroc exposure parameters, C avg 
(AUC divided by dose interval: AUC/t), C min , and ECC, were 
obtained using a Bayesian feedback procedure, as described 
previously. 21 ECC was defined as the constant concentration 
that gave the same average inhibition as the time varying 
concentration over a dosing interval and was calculated 
using an in vivo half-maximal inhibitory concentration (IC 50 ) of 
7.65ng/ml. 15 ECC is theoretically the most appropriate expo- 
sure variable when comparing different dosing frequencies 
(q.d. vs. b.i.d.) and/or different-shaped maraviroc concentra- 
tion-time curves, such as seen with inducers, inhibitors, or 
neutral agents. 23 Exposure variables were calculated exclud- 
ing BLQ concentrations. The BLQ information was reported 
as a separate variable as it was considered to be the result of 
poor adherence rather than insufficient dosing. 

Exploratory analysis. An exploratory graphical analysis was 
performed of the efficacy-concentration data in maraviroc- 
treated patients. The success rate per quartile of exposure 
was also tabulated and plotted. 
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Modeling. Viral RNA, with treatment success defined as VL <50 
copies/ml at 48 weeks, was analyzed as a binary variable using 
GAMs in which the effect of prognostic factors was examined. 

A commonly used statistical model for binary data is the 
logistic regression model. In the frequently used linear logis- 
tic regression model, it is assumed that the log-odds of the 
probability of an event are linear given the prognostic factors 



x i1? x i2 , ...x ip : 



log 



p(yil*ii,*i2>--*i P ) 
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Although attractively simple, these traditional linear logis- 
tic regression models might fail, because real-life effects are 
generally not linear. For this reason, the more general method 
of GAMs 24 has been used for analyzing clinical endpoints in 
HIV patients and their relationship with maraviroc exposure 
and prognostic factors. 

In the generalized additive logistic model, a more general 
functional form can replace each linear term: 

p(y, |x l1 .x i9 ,....x in ) 

;t) -a+<w+-w b 

where f can be an unspecified ("nonparametric") function, 
which is estimated in a flexible manner using an algorithm 
whose basic building block is a scatter plot smoother (e.g., 
locally weighted scatterplot smoothing). A smoother is a tool for 
summarizing the trend of a response variable as a function of 
one or more predictors. It produces an estimate of the trend that 
is less variable than the response variable (y). Smoothing takes 
place by local averaging of y-values of observations having pre- 
dictor values close to the target values. A simple example of a 
smoother is a running mean or moving average. However, in a 
GAM not all of the functions need to be nonlinear, and linear 
parametric forms can be easily mixed with nonlinear terms. 

The GAMs were built using the automated step-wise search 
developed in S-PLUS (version 6.2; Insightful, Seattle, WA). 
The automated step-wise search selects the best GAM using 
alternating forward selection and backwards deletion given the 
range of models. A series of candidate relationships (e.g. linear, 
log-transformation, spline, locally weighted scatterplot smooth- 
ing) that describe how each particular prognostic factor might 
enter the model is defined for every prognostic factor. The final 
models were built by evaluating all the candidate forms for each 
prognostic factor in a step-wise manner (Supplementary Data 
online has Splus code for the final-step GAM model and a spec- 
imen part data set). The prognostic factors considered for inclu- 
sion in the GAMs are listed in Supplementary Table S1 online. 

The uncertainty in the GAMs for efficacy endpoints was 
quantified using a bootstrap technique. 25 The data set for 
the GAM analysis of the clinical efficacy endpoint was boot- 
strapped 1 ,000 times. For every bootstrap sample, the relevant 
final GAM was fitted providing 1 ,000 sets of GAM parameters 
for each clinical endpoint and the median; 95% confidence 
intervals of the predicted clinical endpoints for the 1 ,000 boot- 
strapped data sets were compared against the observed clini- 
cal endpoints. The predictive performance of the final GAMs 
was judged acceptable if the observed values fell within the 
predicted distributions. 



A description of how trial simulations were performed is 
available in Supplementary Data online (Trial outcome 
simulations). 
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Study Highlights 

WHAT IS THE CURRENT KNOWLEDGE ON THE 
TOPIC? 

y The outcome of long-term HIV-1 treatment may 
be influenced by a number of disease and treat- 
ment factors as well as exposure to test treat- 
ment. Generalized additive models (GAM) have 
proved useful in assessing multiple linear and 
nonlinear relationships in landmark analyses in 
human immunodeficiency virus (HIV) and hepa- 
titis C virus studies. 

WHAT QUESTION DID THIS STUDY ADDRESS? 

y This analysis assessed exposure-response 
relationship for maraviroc (+ optimized back- 
ground therapy) at 48 weeks in the context of 
other predictive factors (including adherence) in 
heavily pretreated HIV-1 -infected patients. 

WHAT THIS STUDY ADDS TO OUR KNOWLEDGE 

y No maraviroc concentration-based exposure- 
response relationship was identified once other 
prognostic factors, including categorical vari- 
ables of treatment group (maraviroc vs. placebo) 
and adherence (drug holiday), were included in 
the model as strong predictors. 

HOW THIS MIGHT CHANGE CLINICAL 
PHARMACOLOGY AND THERAPEUTICS 

S Although GAM analysis is an efficient method of 
assessing covariates in large datasets, it is also 
important to carefully consider how adherence 
is factored into population pharmacokinetics 
with sparse sampling and how poor adherence 
is used in exposure-response relationship and 
subsequently in clinical interpretation. 
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